arXiv: 1504.00676v2 [astro-ph.GA] 22 Sep 2015 


Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 24 September 2015 (MN style file v2.2) 


Direct collapse black hole formation via high-velocity 
collisions of protogalaxies 


Kohei Inayoshi^* f, Eli Visbal^J§, Kazumi Kashiyama^^H 

^Department of Astronomy, Columbia University, 550 W. 120th Street, New York, NY 10027, USA 

^Department of Astronomy; Theoretical Astrophysics Center; University of California, Berkeley, Berkeley, CA 9^720, USA 


24 September 2015 


ABSTRACT 

We propose high-velocity collisions of protogalaxies as a new pathway to form su- 
permassive stars (SMSs) with masses of ^ 10® Mq at high redshift {z > 10). When 
protogalaxies hosted by dark matter halos with a virial temperature of ^ 10"^ K col¬ 
lide with a relative velocity > 200 km s“^, the gas is shock-heated to ~ 10® K and 
subsequently cools isobarically via free-free emission and He"*", He, and H line emis¬ 
sion. Since the gas density (> 10^ cm“®) is high enough to destroy H 2 molecules by 
collisional dissociation, the shocked gas never cools below ^ 10^ K. Once a gas cloud 
of ~ 10® Mq reaches this temperature, it becomes gravitationally unstable and forms 
a SMS which will rapidly collapse into a super massive black hole (SMBH) via general 
relativistic instability. We perform a simple analytic estimate of the number density of 
direct-collapse black holes (DCBHs) formed through this scenario (calibrated with cos¬ 
mological N-body simulations) and find udobh 10“® Mpc“® (comoving) by z = 10. 
This could potentially explain the abundance of bright high-z quasars. 

Key words: black hole physics, cosmology: theory, cosmology: dark ages, reioniza¬ 
tion, first stars, galaxies: formation, quasars: supermassive black holes 


1 INTRODUCTION 

The e xistence of br i ght q u asars (QSOs) at z > 6 — 7 IIFarJ 
2006; IWillott et al.l I2OI0I : iMortlock et al.l I2OIII : IWu et ahl 
20151 ) presents an intriguing question: how do supermas¬ 
sive black holes (SMBHs) with masses > a few x 10® M© 
form within the first billion years after the Big Bang? 
Perhaps the simplest possible explanation is that the 
10 — 100 M© black hole remnants of the first genera¬ 
tion of stars grow into these supermassive black holes via 
gas accretion. However, this requires essentially uninter¬ 
rupted Eddington limited accretion for the entire history 
of the Universe, which seems unlikel y due to radiative 


feedback from the accreting BH (e.t 

r. I.Tohnson Kr. Bromm 

2007; 

Milosavlievic. Couch & Bromml 

2 OO 9 I: Park & Ricotti 

2011, 

2 OI 2 I: Tanaka. Perna & Haiman 

2012tl. Maior merg- 


ers of BHs are not expected to accelerate growth signif¬ 
icantly. This is because the kick velocity of the merged 
BH is typi cally larger than the escape velocity of it s host 
galaxy fe.g. lHerrmann et al.ll200^ : lKoDDitz et al.ll2007l ) lead- 
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ing to ejection and ha lting gas accretion (but see also 
iTanaka fc Haimanll2009l) . 


An alternative scenario is the formatio n of supermas- 
sive s tars (SMSs) with a ma ss of > 10® M 
1994; | BrpmnL_&LLoeb| 20031: Koushia ppas, Bullock fc Dekel 


2004 : iBegelman. Volonteri fc R.eesl 1200^ which directly 

collaps e into SMBH s via general relativistic insta¬ 
bility (IChandrasekharl 19641 : IZeldovich fc Novikov! Il97ll : 
IShibata fc Shanirol 12002 ). The larger mass of these seed 
BHs reduces the accretion time necessary to reach the 
masses implied by high-redshift QSOs. Formation of a 
SMS requires a ~ 10® M© metal-poor gas cloud with 
no molecular hydrogen (H 2 ) in a massi ve dark matter 
halo with virial temperature of ~ 10^ K (iBromm fc LoebI 
I2OO3I) . In the absence of H 2 , the gas cloud can only 
cool by atomic hydrogen (Lya, two-photon, and H~ 
free-bound em issions) and the temperature remains ~ 
8000 K (e.g., lOmukail I2OOII) . Once such a massive gas 
cloud is assembled, it collapses monolithic ally and isother- 
mally without significant fragmentat ion (|Brqnim__fc_j 2 Qeb 
2OO3I : Shang. Bryan fc Haimad 201^; Regan&Haehneh 


2009b|yjLatifet_aH 20131 : Inavoshi, OmuEi fc Tasked 20141 : 
' l2015i). 


Becerra et al.l l2015l ). After the collapse, a single proto¬ 
star with a mass of ~ 1 M© is formed at the cen¬ 
ter of the cloud. The p rotostar grows via rapid gas ac- 
cretion, > 1 M© yr“^ (llnavoshi. Omukai fc Tasker l2014l : 
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iBecerra et al.l 2015ll. and become s a SMS within ^ 1 
Myr jHosokawa, Omu kai fc Yorkd I2OI2I : iHosokawa et alJ 
l2013l : ISchleicher et al.ll2013l 'l. 

Suppressing molecular hydrogen cooling is the largest 
obstacle to high-redshift SMS formation. One way to 
achieve this is through H2 photodissociation by far- 
ultraviolet (FUV) photons in the Lyman-Werner (LW) 


2002 

: iBromm Sz LoebI 12003 

: IShane:. Brvan Sz HaimanI 

201c 

: Inavoshi fc Omukail 

2011: 

Latif et al.l 20141: 

Sueimura. Omukai fc Inouel l2014l!. In 

order to dissoci- 


ate H2 in a massive dark matter halo, the required FUV 
intensity is ~ 1500 (in units of 1 0~^^ erg s“^ cm~^ 
sr“^ Hz~^: ISuaimura, Omukai fc Inoua l2014l 'l. Given that 
star-forming galaxies are also st rong X-ray sources , Jr,w 
may increase by a factor ~ 10 Jlnavoshi fc Omukail I2OI1I : 
llnavoshi fc Tanakal |2014| '| . This is due to the electron- 
catalyzed reactions (H -f e“ —>■ H“ -|- 7 ; H~ -f H ^ H2 -I- e“), 
which promote H2 formation in low-metallicity gas. Massive 
dark matter halos irradiated with ~ 10"^ are expected 
to be extremely rare. Direct-collapse BHs (DCBHs) formed 
in this way may not be able to account for the abundance 
of ob s erved high-z quasa r s llPiikstra. Ferrara fc Mesingerl 
I2OI4I : llnavoshi fc Tanakal 12014 1. We point out that a 
variation of this scenario based on the synchronized 
formation and merger of pairs of Tvir ~ 10^ K dark 
matter halos, may produce enough DCBHs to explain the 
observations dVisbah Haiman fc BrvanI l2014bf l. However 
additional numerical simulations are required to validate 
this possibility. 

A different pathway to H2 suppression is collisional dis¬ 
sociation (H2 + H —>■ 3H), which can occur if the metal 
poor gas reaches high density and temperature, satisfying 
(n/10^ c m~^) X (T/IO^ K) ^ 1 (th e so-called “zone of no 
return”!. llnavoshi fc Omukail (l2012l ! proposed that galactic- 
scale shocks can satisfy this condition. If the shock happens 
at the central region of a massive halo <0.1 f?vir, the density 
and temperature of the shocked gas become n > 10"^ cm“^ 
and T > 10^ K, and efficient collisiona l dissociation of 
H2 ca n occur. However, the simulations of iFernandez et ^ 
(I2OI4I! showed that for several examples of less massive halos 
with Tvir < 10^ K, shocks do not reach the center preventing 
SMS formation. It was also pointed out that in typical ha¬ 
los (not the high-velocity collisions discussed in this paper), 
the zone of no return cannot be reached without radiative 
cooling, w hich may lead to star formation a nd prevent SMS 
formation (IVisbal. Haiman fc Brvanll2014ah . SMS formation 
may still be possible in larger halos (Tvir > 10^ K) if shocks 
reach the centers of the halos before significant amounts of 
stars of formed. This requires further study with numerical 
simulations. 

Another proposed SMS formation scenari o is b ased 
on massive-galaxy mergers dMaver et al.l I2OI0I . l2014h . A 
merger can drive strong gas inflow and supersonic turbu¬ 
lence in the inner galactic core, which prevents significant 
fragmentation of the gas even wit h some metals (but see 
iFerrara. Haardt fc Salvaterrall20l(j !. If efficient angular mo¬ 
mentum transport can be sustained in the inner ~ 0.1 pc for 
a sufficiently long time (which requires confirmation from 
further numerical simulations) a SMS of > 10® Mq could 
form. DCBHs from such SMSs might explain the observed 
abundance of high-2 QSOs. 


In this paper, we propose high-velocity collisions of pro¬ 
togalaxies as a new pathway to form SMSs and DCBHs 
at high redshift. As observed in the local Universe, a frac¬ 
tion of galaxies and also clusters of galaxies collide with 
a much larger velocity than the typical pe culiar velocity 
(e.g.. the Taffy galaxy and the bullet c l uster: ICondon et al.l 
]T9^_ CondrH]j _He2ou^ Jarretd I2OO2I : I Tucker et al. 19981 : 


Markevitch et al.l 2002l! . At the interface of such collid¬ 


ing g alaxies, shock-induced s t arbursts have been confirmed 
fe.g.. iLarson fc Tinslevl[l97l : ISaitoh et al.l|2009l! . We show 
that, when a similar collision with a high-velocity > 200 km 
s“^ happens between metal-poor galaxies, a hot gas (~ 10® 
K) forms in the post-shock region and the subsequent ra¬ 
diative cooling makes the gas dense enough that any H2 
molecules are destroyed by collisional dissociation. Once a 
gas clump of ~ 10 ® Mq with a low concentration of H2 due to 
collisional dissociation forms, its gravitational collapse can 
be triggered and a SMS forms. Note that our scenario does 
not require supersonic turbulence or extremely efficient an¬ 
gular momentum transfer as in the galaxy merger scenario. 
We estimate the abundance of SMSs and DCBHs produced 
by high-velocity galaxy collisions, and show that it can be 
comparable to that of high-2 QSOs. 

This paper is organized as follows. In §2, we derive the 
necessary conditions to form SMSs in protogalaxy collisions. 
In §3, we estimate the number density of protogalaxy col¬ 
lisions resulting in the SMS formation, and show that the 
DCBHs from such SMSs could be the seeds of high-redshift 
QSOs. Finally, we summarize and discuss our results in 
§4. Throughout we assume a ACDM cos mology consisten t 
with the latest constraints from Planck jAde et al.ll2014l !: 
Da = 0.68, Dm = 0.32, Db = 0.049, h = 0.67, erg = 0.83, 
and ris = 0.96. 


2 SMS FORMATION VIA PROTOGALAXY 
COLLISIONS 

Generally speaking, a SMS can form from a > 10® Mq 
metal-poor gas clump without H2. Since H2 can form via the 
electron-catalyzed reactions (H -|- e“ —^ H^ -|- 7; H~ TH—^ 
H2 -I- e“), it must be efficiently dissociated. We propose a 
SMS formation scenario where H2 is dissociated from the 
shocks produced by a high-velocity collision of two dark 
matter halos. In this section, we show that SMS formation 
requires the relative velocity of the colliding protogalaxies 
to be in a specific range. If the collision velocity is too low, 
the gas will not be shocked to sufficient temperature and 
density to dissociate H2. On the other hand, if the velocity 
is too high and the shock too violent, the gas will be dis¬ 
rupted before it can cool via atomic hydrogen and form a 
SMS. This velocity window depends on redshift due to the 
evolution of the typical gas properties of pre-shocked gas 
within dark matter halos. 


2.1 Protogalaxy properties 

Next, we describe the properties of dark matter halos and 
the gas within them as a function of redshift. This sets 
the collision velocity bounds for SMS formation which are 
derived below. We consider protogalaxies hosted by dark- 
matter halos with virial temperatures of Tvir ~ 10^ K, cor- 
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responding to the atomic cooling threshold. Larger halos 
undergo radiative cooling which triggers star formation and 
metal enrichment, inhibiting SMS star formation (see Sec. 
14.2.11 for further discussion). On the other hand, for halos 
Tvir 10"^ K, the small gas mass is not sufficient to form a 
SMS (see equation [SJ. 

The virial mass of a dark matter halo is given by 

(1, 

and the virial radius by 

«*=.350/.-pc (2) 

dBarkana fc Lo^l200ll L 

Simulations show that before cooling becomes efficient 
the central regions of dark matter halos contain a gas core 
with approximately constant d ensity and radius (see e.g. 
IVisbal. Haiman fc Brvanll2014al ') 

7?core ~ 0.1 f?vir. (3) 

The gas core is surrounded by an envelope with a density 
profile roughly given by oc r~^. The entropy profile, defined 
as K — , also has a core with K/K^it 0.1. Here 

Kvii = and hb is 200 H~^ times mean numbe r 

density of baryons (e.g., IVisbal. Haiman &: Brvaiil l2014al i. 
Since T ~ T^ii in the (pre-shock) gas core, the gas density 
can be estimated as 

«o-22cm-^(^i^) . (4) 



Figure 1. The collision-velocity window for supermassive star 
formation in protogalaxy collisions (shaded light-blue region). 
The solid curve shows the lowest velocity required to induce colli- 
sional dissociation of H 2 (equation^. The dashed curve shows the 
highest velocity that meets the radiative shock condition (equa¬ 
tion m. 


limit of the collision velocity for SMS formation, which is 
given by 

Vo > 210 km s“^ (“ 1 ^) ' 

2.3 Upper velocity limit 


The total core gas mass is 

r n ■ \ / 1 4- r \ 

~ 3.0 X 10® Mo j j , (5) 

which is xxp 10 per cent of the total gas mass inside the dark- 
matter halo. 

2.2 Lower velocity limit 


The upper collision velocity limit leading to SMS forma¬ 
tion is set by the radiative shock condition. If the collision 
velocity is too large, the shocked gas starts to expand adia- 
batically before the radiative cooling sets in. In this case, the 
shocked gas cannot become dense enough to dissociate H 2 
collisionally, and the shock-induced SMS formation cannot 
be triggered. Therefore, the radiative cooling time needs to 
be shorter than the dynamical time of the shock. 

The dynamical time of the shock can be estimated as 


The lower collision velocity limit is set by the requirement 
that H 2 is collisionally dissociated. This will happen if the 
shock is strong enough for the gas to enter the so called 
“zone-of-no-return” fsee llnavoshi fc Omukail|2012l . and Ap¬ 
pendix A). This region in temperature-density space is de¬ 
fined by 

for n < 10^ cm“®, where T and n are the post-shock tem¬ 
perature and density, respectively. For a collision velocity, 
Vo, much larger than the sound speed of the pre-shocked gas 
(xxx 10 km s“^). 


T = ~ 8.5 X 10® K 

Idkb 


( _ ^2 _f (7) 

V250 kms-i/ ’ ^ ’ 


and 


n = 4no, ( 8 ) 

where we set the mean molecular weight as p = 0.6. From 
equations ( 0 , ®, 0 , and 0 , one can obtain the lower 


tdyn — 




Vo 


~ 2.0 X 10® yr 


/ Rcore \ 


(50 pc ) 

V250 km s-iy 


( 10 ) 

while the radiative cooling time of the shocked gas is given 
by 


^cool = “;:rT • 

^■^^rad 


( 11 ) 


Here, Arad(n,T) = n^A(r) is the cooling rate (in units 
of erg s“^ cm“®). The cooling function A(T) consists of 
the contributions from atomic H line emission at T ~ 10'* 
K, atomic He*' and He line emissions at T 10® K, 
and the bremsstrahlung emission oc at T > 10® K 

llSutherland fc Dopitall 19931 : iGlover fc JaPDse3200m . In our 
scenario, the shocked gas temperature initially ranges up to 
5 X 10® K < T < 5 X 10® K, which corresponds to the colli¬ 
sion velocity of 150 km s“* S, vo S, 500 km s“* (see equation 
0 . For such gas, the cooling function can be well approxi¬ 
mated by Ao — 5 X 10“^“* erg s“* cm®. Given that n oc T“* 
during the radiative contraction, the cooling time becomes 
shorter for a lower-temperature, i.e. the gas contracts in a 
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thermally unstable way. Thus, the radiative cooling time is 
essentially given by substituting Aq into equation (HU; 


^cool,] 


2.8 X 10"* yr 


/ no / no \-i 

V250 kms-i/ VlOcm-3/ 


( 12 ) 

From equations P-ll, Cq}, and m, the radiative shock 
condition (tdyn ^ icooi) can be rewritten as 


Vo < 620 km s ^ 


Tvir y+zy/" 

104 K j V 16 j 


(13) 


Fig. [T] shows the collision velocity window for SMS for¬ 
mation as a function of redshift (the shaded region). For col¬ 
lisions of Tvir ~ 10"* K halos in this window, gas is shocked 
into the zone-of-no-return and H 2 molecules are destroyed 
by collisional dissociation. Additionally, the shocked gas ra- 
diatively cools to ~ 10"* K within the shock dynamical time. 
Once the > 10® Mq gas cloud is assembled and cools, it 
becomes u nstable due to self-gravity, and a SMS can be 
form ed ie.g. llnavoshi. Qmukai fc Taskeil2014l : lBecerra et al.l 
I2OI5II . Note that since the total mass of the colliding gas in 
dark matter halos with Tvir — lO'^ K can be as large as 
~ 10® Mq, several SMSs may be formed at once in the col¬ 
lisions we consider. These SMSs would result in DCBHs of 
> 10® Mq at z > 10 , and as we show in the following section 
could potentially be the seeds of high -2 QSOs. 


3 DCBH ABUNDANCE FROM 
HIGH-VELOCITY COLLISIONS 

Precisely estimating the number of collisions which result in 
SMSs and DCBHs is very challenging because it depends on 
detailed nonlinear physics. This most likely necessitates N- 
body simulations, however, the rarity of these collisions (we 
estimate ~ 10~® Mpc“® from 2 = 10 — 20) requires simula¬ 
tions much larger than are feasible with current computers. 
Here we address this issue by performing a simple order of 
magnitude estimate with an analytic formula based on ideal¬ 
ized assumptions and calibrated with a N-body simulation. 
We find that the number density of DCBHs could be high 
enough to explain observations of high -2 QSOs. However, we 
emphasize that our estimate has large uncertainties which 
we discuss in §4. 



Figure 2. The peculiar velocity PDF, p{v), measured from our 
N-body simulation (solid curves) and the best fits described in 
Section 3 (dashed curves) tor redshifts 2 = 20, 15, 12, 10 (from 
left to right). 


threshold, Ufast is the velocity of the fast-moving halo neces¬ 
sary to form one (or several) SMS(s), rifast = /fastUh is the 
number density of halos with peculiar velocity greater than 
this value (but below the maximum value), and h is the im¬ 
pact parameter required for SMS formation. Note that these 
values are all initially calculated in physical units and the 
abundance is later converted to comoving units to compare 
with QSO observations. We compute the halo number den¬ 
sity w ith the Sheth-Tormen mass function (ISheth fc Tormeiil 
I 1999 I) and consider a mass range of (0.5 — 1) x Mcooi (equa¬ 
tion [T]). For the impact parameter, b, we use ten per cent of 
the virial radius (equation|3}- We assume Ufast = Umin —Utyp, 
where utyp is the typical halo peculiar velocity (assumed to 
be 40 km s“^) and Wmin is shown in Fig. [T] (solid curve). The 
fraction of fast-moving halos is given by 

r'^max — ^typ 

/fast = / p{v)dv, (15) 

where Umax is given in Fig. [T] (dashed curve) and p{v) is the 
peculiar velocity probability density function (PDF) which 
we estimate from a N-body simulation as described in the 
following subsection. 


3.1 Collision rate 

We estimate the high-velocity protogalaxy collision rate by 
considering the number of dark matter halos just below the 
atomic cooling threshold that collide with a relative velocity 
in the range shown in Fig. [T] For simplicity, we consider 
one halo moving with a very high peculiar velocity and the 
other with a typical velocity (determined with the N-body 
simulation described below) in the opposite direction. 

Making the idealized assumption that halo positions 
and velocities are randomly distributed (i.e. ignoring clus¬ 
tering and coherent velocities, which are discussed in Sec. 
14.2.31 . the collision rate per volume is given by 

dUcoll ,2 j. 2 i_2 

— ^fast^h ^ '^fast ^ — /fast^h ^ “^fast ^ 

where rih is the number density of all halos near the cooling 


3.2 N-body simulation and velocity PDF 

To estimate the dark matter halo peculiar velocity PDF, 
we ran a cosmological N-bod y simulation w ith the pub¬ 
licly available code GADGEt2 llSpringelll2005h . The simula¬ 
tion has a box length of length 10 h~^ Mpc (comoving) and 
mass resolution of 768® particles, corresponding to a par¬ 
ticle mass of 1.96 x 10® /i“® Mq. Snapshots were saved at 
2 = 20, 15, 12, and 10. We used the rockstar halo hnder 
llBehroozi. Wechsler fc ^i^l2013l ) to locate dark matter ha¬ 
los and determine their masses and velocities. 

In Fig. [21 we plot the velocity PDF of M = (0.5 — 
1) X Mcooi dark matter halos at 2 = 20, 15, 12, and 1 0. We 
also plot fits with a f orm g uided bv ISheth fc Diaferid (1200 ll) 
and iHamana et al.l (l2003l ) who argue that the PDF is a 
Gaussian distribution in each velocity component (Maxwell- 
Boltzmann distribution for the total ID velocity) at fixed 
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Figure 3. The estimated number density of DCBHs created t\ 
cumulative number density as a function of redshift (right panel). 



redshift z 

high-velocity collisions per unit redshift (left panel) and the 


halo mass and local overdensity. For halos near the cooling 
threshold this leads to 


_ f dSp(S)(l + bhS}p(vlS) 
J dSp{S){l + bi,5) 


(16) 


where p{S) is the cosmological overden sity PDF, bh is 
the S heth-Tormen dark matter halo bias (ISheth fc TormenI 
I1999I I. and p(w|(5) is the velocity PDF at fixed <5. The over- 
den sity PDF is assumed to be a lognormal distribution (see 
e.g. I Coles fc Jonedfig^ 


p{S) = 


1 

\/27l-cr| 


exp 


[ln(l + J)+a|/2]^ \ 1 

2cr| J 1 + S 


(17) 


The velocity PDF at fixed overdensity is assumed to be a 
Maxwell-Boltzmann distribution with variance 


a^ = {l + Sf'^al (18) 

We set cr| = ln[l + 0.25/(l + 2)] llHamana et al ] |2003l l (which 
determines the size of the region corresponding to S) and fit 
two parameters to our data, a„ and p. For redshifts of 2 = 
20, 15, 12, and 10, we find cr„ =24.04, 27.27, 30.91, and 33.54 
and p = 0.8687, 1.2404, 1.0949, and 1.2081, respectively. 


3.3 DCBH number density 

Using the best fit p{v) discussed above, we find the num¬ 
ber of DCBHs produced from 2 = 10 — 20. We assume that 
the LW background at 2 < 20 suppresses star formation 
and subsequent metal enrichment in halos below the atomic 
cooling threshold. At higher redshift we assume that star for¬ 
mation in minihalo progenitors prevents DCBH formation. 
The total number density formed as a function of redshift is 
given by 

nDCBH(2)=|d2^^^. (19) 

To get the velocity PDF at intermediate redshifts between 
our simulation snapshots, we linearly interpolate the results 
from Fig. [2] In Fig.|31 we plot the number density of DCBHs 
as a function of redshift. Most DCBHs come from high red¬ 
shift due to the lower minimum velocity given in Fig.[T] We 


find a total density by z = 10 of 10“® Mpc“^ (comoving). 
Thus, it seems possible that these DCBHs could potentially 
explain the abundance of high-2 QSOs. 


4 DISCUSSION AND CONCLUSIONS 

We have shown that high-velocity collisions of metal-poor 
galaxies may result in the formation of supermassive stars 
(SMSs). When dark matter halos with a virial temperature 
~ 10'* K collide with a relative velocity > 200 km s“*, gas is 
heated to very high temperature (~ 10® K) in the shocked 
region. The shocked gas cools isobarically via free-free emis¬ 
sion and forms a dense sheet (> 10"* cm“®). In this dense 
gas, H2 molecules are collisionally dissociated, and the gas 
never cools below ~ 10"* K. Such a clump of gas with mass 
~ 10® Mq, once assembled, becomes gravitationally unsta¬ 
ble and forms a SMS, which would directly collapse into a 
black hole (DCBH) via general relativistic instability. We es¬ 
timated the abundance of DCBHs produced by this scenario 
with a simple analytical argument calibrated with cosmo¬ 
logical N-body simulations and found a number density of 
~ 10“® Mpc“® (comoving) by 2 = 10. This is large enough 
to explain the abundance of high-redshift bright QSOs. 


4.1 Observational Signatures 

Next, we briefly discuss the possible observational signa¬ 
tures of SMSs formed through high-velocity collisions of 
protogalaxies. The temperature of the shocked gas in the 
collisions we consider above is T < 10® K (equation [7]). 
The gas cools initially via bremsstrahlung, then atomic He*' 
and He line emissions, and finally atomic H line emission. 
The intrinsic bolometric luminosity can be estimated as 
Lboi ~ 10Mgas,coreU0^/fcooi,max < 10*® erg s“* for our rep¬ 
resentative case (see equations [5] and uni. Given that the 
colliding gas in dark matter halos is mostly neutral, the 
cooling radiation is reprocessed into various recombination 
lines, e.g., Lya, Ha, and He II A1640. The Ha and He H 
A1640 emission lines are particularly interesting, since the 
intergalactic medium would be optically thin to them. If 
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~ a few per cent of the bolometric l uminosity goes into th ese 
lines, which is reasonable (see, e.g.. Ijohnson et ^l201ll . for 
numerical simulations of cooling radiation from hot metal- 
poor gas with ~ 10® K), the emission could be detected 
from 2 < 15 by the Near-Infrared Spectrograh (NIRSpec) 
on the James Webb Space Telescope (JWST) with an expo¬ 
sure time of > 100 h. Due to the high cooling temperature, 
the ratio of He II A1640 to Hq flux is expected to be large, 
which may make it distinguishable from other sources (e.g. 
population III galaxies). In principle it may also be possible 
to detect H Lya emission from protogalaxy collisions. Lya 
emission could constitute a large fraction of the bolometric 
luminosity (perhaps 10 per cent). However, even if a col¬ 
lision is observed in Lya it may be difficult to distinguish 
from other objects such as accreting massive dark matter 
halos. 

A detailed radiative transfer calculation is required to 
accurately predict the emission spectrum, which is beyond 
the scope of this paper. Even though protogalaxy collisions 
may be bright enough to observe, recent collisions are ex¬ 
pected to be extremely rare. At most there will be ~ a few 
in the whole sky, given the event rate, dnco\\/dt{z = 15) ~ 
10~^^ Mpc“® (comoving) Myr“^ and the emission duration, 
~ 0.1 Myr. Thus, it will be extremely challenging to detect 
the signal described above in the near future. 

4.2 Impact of assumptions 

Here we discuss some of the key assumptions we made, and 
how changes to these assumptions would affect our results. 


4-2.1 Metal enrichment 


In §2, we calculated the thermodynamics of the shocked 
gas after protogalaxy collisions assuming zero metallicity. 
This assum ption is valid for gas me tallicity smaller than 
< 10“® Zq lllnayoshi &: Omukaill201^ 1. If the metallicity is 
higher than this critical value, the shocked gas can cool down 
to below ~ 10"* K via metal-line emissions (Cn and Oi) and 
fragment into clumps of ~ 10 Mq, preventing SMS forma¬ 
tion. 


In general, the metal-enrichment of gas in massive dark 
matter halos proceeds in two different ways. The first is in¬ 
ternal enrichment by in situ star formation. Although not 
yet completely understood, the earliest star formation is ex¬ 
pected to be triggered by H2 cooling in progenitor “mini¬ 
halos” with Tvir < K, which eventually assemble into 
the more massive dark matter halos we consider in this pa¬ 
per. The level of self enrichment in minihalos is sensitive 
to t he initial mass function (IMF) of population HI star s 

fe.g. lHirano et ^l2014l : ISusa. Haseeawa fc Tominagall20T3l . 

If the IMF is extremely top heavy, the metal enrichment 
is predominately provided by pair-instability SNe. In this 
case, the metallicity at the gas core inside the dark matter 
halo could be as large as 10~^ — 10“® Zq at « ~ 10 (e.g. 
iGreif et al.ll201ol : IWise et al.l2012|l . On the other hand, if the 
IMF is mildly top heavy, core-collap se SNe from ~ 40 Mp 
stars would be the dominant source jH osokaw a et al.ll201ll : 
IStacv. Greif fc Brommir2012l : iHosokawa et al.l 20121 1. In this 
case, the metallicity may be one o rder of magnitude lower 
(5, 10~^ Zq) at the same redshift jHeger fc VWosl^l2002l : 
iNomoto et al.lli^OOd l. 


In our abundance estimates of SMS formation through 
high-velocity collisions, we only consider dark matter ha¬ 
los with Tvir < 10“^ K. We assume that below 2 = 20 
the abundance of H2 required for star formation in mini¬ 
halos i£_sufficiently_su£pr^sed_W LW background radiation 
fe.g. iHa iman. Abel fc Reesll200ol: iMachacek^ Br van fc Abell 


I2001I : Wise fc Abel 20071 : O’Shea fc Norman 2008ll . The re- 
quired LW background flux for this to occur is estimated 
to be J lw ^ 0.2 — 2 (3 x 10~ ^ — 4 x 10“^) at 2 = 
15 (20) (IVisbal. Haiman fc Brvanll2014bll . The anticipated 
LW background flux is Jlw ~ 0.1 — 10 (0.01 — 20) at 2 ~ 
15 (201 fe.g. lA hn et al. l2009l : [Johnson. Dalla fc Khochfail 
l20ld : IVisbal et al.l 2014 1. which depends on the detailed 
properties of Pop HI stars and the efficiency with which they 
are produced. While there are certainly large uncertainties 
in the LW background, we find our assumption of minihalo 
star formation suppression to be reasonable. 

The second way in which halos can obtain metals is 
through external enrichment by galactic winds from nearby 
massive galaxies. Semi-analytic models predict that the in- 
tergalactic medium can be polluted by this effect lead¬ 
ing to a n average metallicity of (Z) ~ 10~^ Zq by 2 > 
12 (e .g. iTornatore. Ferrara fc Schneideil l2007l : iMaio et al.l 
I2OIOII . However,the fraction of the intergalactic medium that 
has been polluted is expected to be small at the high red- 
shifts important for our calculation (e.g th e estimated vol¬ 
ume filling factor is ~ lO”'^ for 2 > 12; IPallottini et al.l 
l2014h . Thus, external metal enrichment is unlikely to im¬ 
pact our assumption of zero metallicity. 

In summary, it is reasonable to neglect the effect of 
metal cooling for dark-matter halos with a virial temper¬ 
ature Tvir JS 10"* K at 2 > 10 after the LW background 
suppresses star formation in minihalos. 


4-2.2 Gas thermodynamics 

In this paper, we derived the conditions for SMS formation 
in protogalaxy collisions (Fig. [T]) based on the “zone of no 
return” shown in Appendix A. This zone is obtained from a 
one-zone calculation of thermodynamics of the shocked gas. 
Of course, galaxy collisions are actually three-dimensional 
phenomena; detailed hydrodynamical simulations are neces¬ 
sary to confirm our scenario. 

We obtained equation (j6]) by assuming that the shock 
is plane-parallel and steady. This assumption would be valid 
for nearly head-on collisions and timescales shorter than the 
shock dynamical timescale. Accordingly, we set the maxi¬ 
mum impact parameter as 6 ~ 0.1 i?vir, which corresponds 
to the size of the gas core of an atomic-cooling halo. How¬ 
ever, galaxy collisions occur typically with a larger impact 
parameter b ~ Rvir. A critical impact parameter for SMS 
formation needs to be identified by numerical simulation of 
protogalaxy collisions. We note that the formation rate of 
SMSs and DCBHs in our scenario is somewhat sensitive to 
this critical value (oc h^). 

As mentioned in §2, the shocked gas in the zone of no 
return is thermally unstable. Once the instability is trig¬ 
gered, fluctuations in the shocked gas grow and form c lumpy 
structures with a length scale of < Cstcooi (iFie Idl ll965lL As a 
result of this, the structure of the shocked gas deviates from 
the plane-parallel sheet in a cooling time. Unfortunately, 
our one-zone calculation cannot capture these effects. Note 
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that, as long as the H 2 abundance is suppressed, the cool¬ 
ing length is kept shorter than the Jeans length, thus the 
thermal instability does not necessarily result in a smaller 
fragmentation mass. Nevertheless, the effects of the ther¬ 
mal instability on SMS formation need to be studied using 
high-resolution simulations. 

We implicitly assumed that after radiative cooling of the 
shocked gas becomes irrelevant (i.e. fcooi ^ iff), the corre¬ 
sponding Jeans mass is assembled, perhaps within ~ ts, and 
the gas clump collapses due to its self-gravity. Our one-zone 
calculation cannot address how the mass assembly process 
actually proceeds in detail. Even when the mass budget is 
large enough > 10® Mq , the mass assembling may be halted, 
e.g., due to the angular momentum of the gas. A detailed 
numerical simulation is also required to clarify this point. 

Subsequent mass accretion onto the DCBHs formed 
in protogalaxy collisions is also uncertain at this stage. 
This needs to be clarified in order to address whether such 
DCBHs can be the seeds of high-z QSOs. The initial mass of 
the DCBHs is ~ 10® Mq whereas the total gas mass of each 
colliding galaxy is > 10® Mq. We also note that the DCBH 
is unlikely to be hosted by the dark matter halo, at least 
just after the formation, because the collision velocity of the 
parent halos is much larger than the virial velocity. Nev¬ 
ertheless, continuous mass accretion from the intergalactic 
medium may be expected since high-velocity collisions typ¬ 
ically occur in over-dense regions. Additional galactic and 
intergalactic-scale calculations including radiative feedback 
from accreting BHs are required to confirm this. 

4-2.3 High-velocity collision rate estimate 

There are a number of uncertainties associated with the var¬ 
ious assumptions we made to estimate the number density 
of DCBHs produced from high-velocity protogalaxy colli¬ 
sions. First, we note that our estimate depends strongly on 
the precise values of Umin. Due to the steepness of p{v) at 
large v, a 20 per cent decrease in Umin increases the abun¬ 
dance of DCBHs by more than an order of magnitude. The 
abundance also depends strongly on the value of the impact 
parameter needed to create a black hole (nDCBH oc b^). Fu¬ 
ture hydrodynamic simulations of individual collision events 
are needed to constrain Umin and b. 

Additionally, the small box size of our simulation may 
systematically reduce the abundance of halos with high pe¬ 
culiar velocity. This is because large-scale density fluctua¬ 
tions, corresponding to modes larger than the box are ar¬ 
tificially removed. We leave it to future work to determine 
how much this effect could enhance the number density we 
compute here. 

We considered the case of one fast-moving halo and one 
halo at typical peculiar velocity in the opposite direction. Of 
course there can be other combinations of peculiar velocities 
and collision angles which lead to a DCBH. We find that we 
get similar results performing the more complicated analysis 
of adding up the contribution from all angles and different 
combinations of peculiar velocities. There is a factor of a few 
enhancement compared to the simple calculation discussed 
above. 

We note that our idealized assumptions of random posi¬ 
tions and velocity directions are not expected to be accurate 
and estimate their impact on our number density estimate 



Figure 4. The velocity PDF computed directly from our N-body 
simulation at z = 15 (solid curve) and computed after subtract¬ 
ing away the local coherent velocity as described in Section 4.2.3 
(dashed curve). 


here. We expect that fast halos will preferentially be found 
in over dense regions, possibly regions which will soon virial- 
ize. The halo density enhancement (ignored in our estimate) 
is given by (1 -I- bi^S) and the density of DCBHs will depend 
on this value squared. At z = 15, in a region that is about 
to virialize {5 x Sc — 1.686), the density enhancement of 
DCBHs is roughly 100. 

Our assumption of random velocities most likely over¬ 
estimates the number density of DCBHs formed. This is be¬ 
cause on small scales there will be some velocity coherence 
between nearby halos, reducing their relative velocities. To 
estimate the impact of this effect we recompute p{v) from 
our simulation, and for each halo subtract away the mass- 
weighted mean velocity of all other nearby halos in the mass 
range M = (0.5 — l)iV/cooi. We include all halos within the 
typical separation length of these halos, Rc = . In¬ 

creasing the value of Rc by a factor of two does not sig¬ 
nificantly affect our results. If this distance is taken to be 
significantly smaller there are not enough neighbors to com¬ 
pute the coherent velocity. This PDF at 2 = 15 is shown in 
Fig- SI At high V, it is more than an order of magnitude lower 
than p{v) obtained without subtracting coherent velocities. 
The typical peculiar velocity is reduced by 25 km s“^. The 
relative changes at 2 = 20 are similar. We find that these 
two effects (the high-u p{v) and the typical v) lower the 
abundance of DCBHs by approximately two orders of mag¬ 
nitude, which may roughly cancel when combined with the 
correction due to the density enhancement discussed above. 

Despite the large uncertainties described above, the 
high-velocity collision of protogalaxies is an interesting path¬ 
way to form SMSs and DCBHs without extremely strong LW 
radiation and could explain the abundance of high -2 bright 
QSOs. In future work, we plan to perform detailed numer¬ 
ical studies on the gas dynamics of colliding galaxies and 
the event rate of appropriate collisions to determine if these 
events could truly be responsible for the first SMBHs. 
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APPENDIX A: ZONE OF NO RETURN: 
EQUATION (5) 

Following llnavoshi fc Omukail (l2012fl . here we calculate ther¬ 
mal evolution of the shocked gas flow discussed above. 
We assume that the flow is steady and plane-parallel 
dShapiro fc Kand Il98lt) , which is appropriate for almost 
head-on collisions in the shock dynamical time. The con¬ 
servation of mass and momentum between the density po, 
pressure po and flow velocity vg in the pre-shock flow and 
those in the post-shock flow p, p and V give: 

pV = povo, (Al) 


pV^ + p = povo + Po- (A2) 


After crossing a shock front, the gas looses the thermal en¬ 
ergy via radiation following the energy equation. 


de d / 1 Anet 

dt ^ dt\p) p ’ 


(A3) 


where e is the specific internal energy, d/dt the Lagrangian 
time derivative, and Anet the net cooling rate per unit vol¬ 
ume (in units of erg s“^ cm“^). Assuming a strong shock 
(i.e. poVg 3> Po), the initial temperature and density for the 
post-shock flow are given by equations © and (j8} using the 
velocity and density of the pre-shock flow. 

As long as the cooling time of the shocked gas is shorter 
than the free-fall time (tfi = y^3n/32Gp), which is the 
grow th time-scale for gravitational instability (e.g. iLarsonl 
[lUi), we follow the thermal evolution solving the equations 
(EU and (1^ . When the cooling becomes inefficient and 
tcooi exceeds ts, the contraction of the layer halts and dense 
clouds begin to develop inside the post-shock region. Once 
the mass exceeds the Jeans limit, the cloud collapses owing 
to its self-gravity in a runaway fashion following the equa¬ 
tion 


d£ ^ P_ 
dt tff 


(A4) 


We consistently solve the chemical reaction networks among 
primordial species (H, H2, e“, H"*", H“, He, He^ and 

He”*"^). We consider radiative cooling by free-free emission, 
atomic lines (H, He, He+) and H2 lines as well as chemical 
cooling/heating. 

Fig. El] shows the thermal evolutionary tracks of the 
post-shock gas. The solid curve which starts from the open 


(filled) circle corresponds to the case where the initial con¬ 
ditions of the post-shocks gas are n = 25 (100) cm“® and 
T = 10® K. The corresponding density and velocity of the 
pre-shock flow are no = 6.3 (25) cm“® and vg = 270 km s“^, 
respectively. Initially, hydrogen is fully ionized and helium 
is neutral. For both cases, the gas cools down to ~ 10"* K by 
free-free emission and atomic line emission (He^, He, and 
H). In the case of the lower initial density (open circle), the 
gas temperature decreases further to ~ 300 K by H2-line 
cooling. On the other hand, in the case of the higher initial 
density (filled circle), H2 formation is suppressed by the col- 
lisional dissociation (H -|- H2 3H). The cooling becomes 
irrelevant when tcooi ^ iff. Then, the gas cloud collapses by 
the self gravity once the corresponding Jeans mass is as¬ 
sembled. In the case of the lower initial density, tcooi ^ tg 
occurs at n ~ 2 X 10® cm“® and T ~ 100 K, and the cor¬ 
responding Jeans mass is a few ~ 100 M©. On the other 
hand, in the case of larger initial density, tcooi ^ iff occurs 
at n ~ 4 X lO"* cm“® and T ~ 8000 K, and the correspond¬ 
ing Jeans mass is > Mj ~ 10® M©. Such a massive cloud 
collapses almost isothermally, mediated by H atomic cool¬ 
ing, and forms a proto-SMS at the center without a ma- 
j or episode of fragmentation and subsequent H2 form ation 
lllnavoshi. Omukai fc Taskeill201'i : iBecerra et al.ll^lSl l. 

The shaded region in Fig I A1 1 represents the “zone of no 
return”. If the gas jumps into this region by a strong shock, 
H2 formation is suppressed by the collisional dissociation 
and massive clouds with > 10® Mq are form. The boundary 
is identified numerically (cross symbols), and the dashed line 
is the fitting, T > 5.2 x 10® K (n/10^ cm“®)“^ for n < 
10"* cm“® (equation [6|). 

Finally, we note the dependence of the density range 
of the zone of no return on the reaction rate coefficient 
of H2 collisional dissociation. We adopt the rate given by 
[Martin. Schwarz fc Mand^ (Il996ll . The value of this coeffi¬ 
cient is ~ 30 times higher than that used in some previous 
studies at n = 10® — 10"^ cm~® (e.g . IShaniro fc Kang||l987l : 
iMackev. Bromm fc Hernauistl2003 ). This results in an order 
of magnitude decrease in the density required to enter the 
zone of no return and explains why these previous studies 
did not find gas entering it. 
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